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In this paper I study the energy resolved supercurrent of a junction consisting of a dirty normal 
metal between two superconductors. I also consider a cross geometry with two additional arms 
connecting the above mentioned junction with two normal reservoirs at equal and opposite voltages. 
The dependence of the supercurrent between the two superconductors on the applied voltages is 
studied. 

PACS numbers: 74.80.Fp, 74.50.+r 



The proximity effect between a normal metal and a 
superconductor has been discussed long time ago . Of- 
ten it is simply described by a spatial dependent pairing 
correlation function ^ which decays from the supercon- 
ductor to the normal metal. However, this description 
is too crude to provide a proper understanding of the 
phenomena observed at low temperatures in the meso- 
scopic systems which can nowadays be prepared in the 
laboratories. For example the detailed description of the 
energy dependence of the effective barrier conductance 
and diffusion coefficient is crucial in understanding 
the behavior of the observed conductance between a nor- 
mal metal (N) and a superconductor (S) at low voltages 
and temperatures. [|| ||| 

In this paper we study the spectral current density fig] ] 
(see also of a quasi-one dimensional SNS junction 
in the dirty limit. This quantity (or, more precisely, the 
angular average of the one defined in [^3|) is defined as, 
at energy e and position x, 

Nj{e,x)=<p x N(p,e,x)> (1) 

where N{p, e, x) is the density of states for momentum di- 
rection p at energy e and position x. The angular brack- 
ets denote angular average. This quantity is thus the 
density of states weighted by a factor proportional to the 
current that each state carries (in a certain direction, 
here x), and thus may also be appropriately referred to 
as the current-carrying density of states. This is obvi- 
ously a useful quantity. For example at equilibrium, the 
(number) supercurrent J s can be written as 



J s = -2v f J -Nj(e,x)h (e) (2) 

where /io(e) = tanh^y and Vf is the fermi velocity. The 
factor of 2 includes the contribution from the two spin 
directions. One convenient way to interprete this formula 
JTJIQ (see also (H[@) is to rewrite h = (1- 2n) where 
n(e), the occupation number, is given by the Fermi func- 
tion at equilibrium. For example at T — eqn (|^) can 
be re-written as (using the symmetry Nj(e) = —Nj(—e) 
) 

J s = 2v f I deNj(e) (3) 

J — oo 

and thus can be interpreted as the current due to the 
occupation of negative energy states. This can also be 
regarded as the diamagnetic response of the supercon- 
ductor if one considers the T = state as one containing 
no quasiparticlcs. Similarly at finite temperature 

J S (T) = J s (T = 0) + 

/OO 
deNj(e){n(e,T) - n(e,T = 0)) (4) 
-oo 

and can be interpreted as the sum of the diamagnetic 
current and the correction due to the thermal redistribu- 
tion of quasiparticles. In particular an important source 
of the decrease of the supercurrent as the temperature 
increases is due to the thermal excitations of quasipar- 
ticles from e < to e > states, which carry opposite 
current. 

In the dirty limit, on which this paper will concentrate, 
Nj can be obtained from (see Appendix for details) 
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Nj(e,x) = —fQ(e,x) 
o 



(5) 



where Nt is the density of states in the normal state, I 
is the mean free path, and Q is given by 
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Tr[r 3 (g H dg H - g A dg A )} 



(6) 



Here g R ' A are the angular averaged of the retarded and 
advanced components of the quasiclassical Green's func- 
tion, d represents spatial derivative. The equilibrium 
(number) supercurrent is thus given by 

NfD 



J, = 



deQho(e) 



(7) 



where D = v/l/3 is the diffusion coefficient. For an SNS 
junction with no electron-electron or electron-phonon in- 
teraction in the N region, Q is independent of the position 
x along the junction within that region. 

The behavior of Q is easiest to understand in the limit 
of very short junction ( Ejj = D/L 2 << A, here L is the 
length of the junction and A is the superconducting gap) 
and small phase difference \- I n this case Q should be 
the same as that of a bulk superconductor under a small 
phase gradient. The response of a dirty superconduc- 
tor to a phase gradient or an external vector potential 
is well-known |l7|| . In this case one can show that the 
entire contribution to the supercurrent arises from states 
at e = A, i.e. Q oc 5(e — A). In contrast the ordinary 



density of states is given by N(e) 



Nt 



Under 



f~Ve 2 -A 2 ' 

a small phase gradient, the gap for quasiparticle excita- 
tions persists and in particular there is no contribution 
to Q for energies within this gap. 

An energy gap e g (< A) exists in general also in an 
SNS junction (except phase difference x = 7r )- This 
gap has been studied before in related situations Q |o| . 
Associated with the existence of this (phase dependent) 
gap is a relatively rapid change of g as a function of 
energy (and phase difference). This has made the nu- 
merical calculation somewhat difficult. For convenience 
I will thus mostly concentrate on results where a small 
pair-breaking term 7 has been included in the self-energy 
(see Appendix). 7 is usually chosen to be 0.05A, though 
occasionally results for 7 = will also be shown for com- 
parison. 

The behavior of Q for a relatively short junction is as 
shown in Fig. [|. At small phase differences Q is large 



only for e near A. If 7 were zero then Q would vanish for 
e below a minigap e g . As the phase difference increases, 
the minigap decreases. Correspondingly the region of en- 
ergy where Q is finite also moves down in energy, though 
it remains large in an energy region up to ~ A. 




FIG. 1. Q (in units of 1/L) for a short junction. 
E D = 1.0A. 7 = 0.05A 



For longer junctions, i.e. L >> y/D/ A or equiva- 
lently Ejj << A, the behavior is somewhat different. 
At a given phase difference, the main region of energy 
where Q is significant is no longer of order A. An exam- 
ple for this evolution as a function of increasing length 
is as shown in Fig |2[ For a given phase difference, the 
energy where Q peaks shifts down in energy relative to 
A as L lengthens. This itself may not be surprising, and 
can be understood by analogy with the behavior of en- 
ergy levels under a change in boundary condition in the 
normal state. 



X = n/4 




0.5 1.0 

El A 



FIG. 2. Q for \ — 7r /4 as a function of decreasing Ed- 
7 = 0.05A 
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The more interesting feature is that a negative dip in 
Q appears at higher energies as the junction lengthens. 
For very long junctions, both the peak and the dip of Q 
move to energies of order (a few tens times) Ejj, with 
almost no features left near A (Fig. ||). This negative 
dip has been speculated to exist recently [El| 
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FIG. 3. Q for a long junction. A = IQQEd- This 
result is for 7 = 

In the above I have assumed that the contacts between 
the normal metal N and the superconducting reservoirs S 
are perfect. If potential barriers exist between the N and 
S regions, then Q decreases in magnitude, with a cor- 
responding decrease in the energy where Q peaks. The 
features discussed above survives for moderate barrier 
resistance Rb between N and S. An example of how Q 
evolves as Rb increases is as shown in Fig. ^. 



From the ideas presented above obviously one can af- 
fect the current flowing between the two superconducting 
reservoirs by changing the occupation of the quasiparticle 
states. Temperature is an obvious candidate. This gives 
the well-known reduction of the supercurrent as a func- 
tion of increasing temperature. An alternative way is to 
create a non-equilibrium situation . Here I shall con- 
sider a steady state situation with the advantage that it is 
easy to analyze. The set-up is shown schematically in Fig 
H Geometries closely related to this has been studied be- 
fore p2|-[2^]. However, these references have concentrated 
on different arrangement of voltages and/or other mea- 
sureable quantities. Here I consider the case where the 
superconductors are at equal voltages, chosen to be zero. 
The normal reservoirs are at equal and opposite voltages 
Vat = ±V, I shall study the dependence of the current 
between the superconducting reservoirs as a function of 
V. 
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FIG. 4. Q for x = tt/4, A = \§E D as a function of 
increasing rt,, the ratio of the barrier resistance Rb to that 
of the normal metal, i.e. rb = Rb(2NfDS/L). Here S is 
the area. 7 = 0.05A. 



FIG. 5. The cross geometry. All 'wires' connecting 
the reservoirs are assumed to be quasi-one-dimensional. 

First we should note that the presence of the side arms 
connected to the normal reservoirs affect the behavior of 
Q via the proximity effect. In order to facilitate later 
discussion, I plotted the quantity Q for this spatial ge- 
ometry for the case A = 10Eu for two phase differences 
in Fig. ^. In this example I have assumed that the arms 
between the normal metal and the superconductor are 
symmetric and of equal length (L x = L y = L in Fig || ) 
and area S. Q is finite only for the x arms connecting the 
superconducting reservoirs, and is constant along them. 
Compared with the case without the side arms (Fig. ||), 
we sec that the behavior of Q is somewhat different in the 
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energy region e < Ed- This is because there is now no 
energy gap for quasiparticle excitations for any position 
within the N region on the cross, even for 7 = 0. How- 
ever, for e > Ed the qualitative behavior of Q is almost 
the same as in the case without the side arms, except 
an overall reduction in magnitude. p8[ In particular the 
sign change of Q remains. 
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FIG. 6. The quantity Q (in units of L x ) for the cross 
geometry of Fig. 5. A = 10E D . 7 = 0.05 A. Results for 
7 = are also shown for comparison. 

Obviously if V — a supercurrent I s only flows be- 
tween the superconducting reservoirs, whereas there is 
no current flowing in or out of the normal reservoirs. At 
V 7^ current is in general finite at any position on 
the two arms. I shall denote the currents as I x and I y . 
Neither I x nor I y are position dependent; moreover, the 
current flowing in and out of the normal reservoirs are 
equal on the one hand and those of the superconducting 
reservoirs equal on the other. One can therefore regard 
the current I x (I y ) as simply flowing between the super- 
conducting (normal) reservoirs. I shall thus continue to 
call I x the supercurrent I s . I shall consider how this I s is 
modulated by the voltage V. All results presented below 
are for T = 0. 

I shall concentrate on an example in the most inter- 
esting regime, where A ~ lOEp. The result for dI s /dV 
at Er> = 0.1 A is as shown in Fig. ^. In this parameter 
range dI s /dV at a voltage V is approximately equal to 
— (NfDS)Q at the corresponding energy e = eV. (c.f. 
the corresponding Q in Fig. ||) Also shown is the value 
of I s at the value V, obtained by adding the integral of 
dI s /dV to the equilibrium value of I s . Note in partic- 
ular that for large V, the supercurrent has actually an 



opposite sign from the equilibrium one, thus producing a 
"7r-junction". ( c.f. 
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FIG. 7. dI a /dV versus V for the cross geometry of 
Fig. 5. En — 0.1A, 7 = 0.05A. Also shown are I s as 
functions of V. Energies (eV) are in units of Ed and 7 S 
is in units of NfEoDSL' 1 . 

To understand dI s /dV, it is necessary to know the be- 
havior of the distribution functions for the quasiparti- 
cles. (see Appendix for the technical details). I shall 
denote these functions on the x- (y-) arms as h,Q^{x) 
(^0,3(2/)) etc. Since we are at T = 0, a small change 
of the voltage at V will affect only the occupation num- 
bers at e = ±eV. In Fig || I have plotted the change of 
the distribution functions <5/io,3 a t a relatively low energy 
when V is increased from below to above eV — e. At the 
S-reservoirs (x = ±L x /2) <5/io.3 = by choice, whereas 
at the normal reservoirs ( y — ±L y /2 ) Sh^ — ±1 and 
Sho = — 1. The behavior of Sho^ is easy to understand 
in this low energy limit, where one can ignore the su- 
perflow (Q), the coupling between the diffusion of the 
two distribution functions (M03 = — M30 are small) and 
where the diffusivity for the particles (oc M33) reduces 
to that of the normal state. Thus (see eq([l^) ) 5h^{y) 
is linear in y and dhs(x) ~ 0. Since there is an energy 
gap at the S-reservoir, the effective diffusivity of the en- 
ergy (oc M o) is suppressed near x — ±L x /2. Thus 5ho 
only has small gradients and hence Sho ~ — 1 everywhere 
except near x = ±L x /2. In the language of the more fa- 
miliar occupation number n(e) = (1 — {ho(e) + h^(e)))/2, 
in this e — > limit 5n(y) is linear in y and thus Sn = 1/2 
at (x,y) = (0,0). 6n(x) is almost constant and w 1/2 
near x w and only changes rapidly to near the S- 
reservoirs. A finite e provides a correction to the above 
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picture as can be also seen from Fig. g[ [|30| The values 
of Sn at e = — eV can be obtained by symmetry since 
n(-e) = (1 + h (e) - h 3 (e))/2. At this energy 5n{y) 
changes from at y = —L y /2 to — 1 at y = L y /2 and 
5n(x) ~ —1/2 near the center of the cross. 
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FIG. 8. The distribution functions at e = 0.24E D 
as functions of x/L or y/L for the cross geometry with 
parameters as in the last figure. \ = 

If 8n(x) were exactly ±1/2 at e = ieF and if one ig- 
nores the fact that Sn is actually x dependent, with eqn 
( |l3| ) (or the equivalence of eqn (^) ) for the current it is 
obvious that dI s /dV will be equal to —NfDQS at the 
corresponding energy. In this case then at large V the 
current I s would be exactly zero. However, the actual 
current consists of both the supercurrent and the con- 
tributions from the gradients of distribution functions, 
(see eq(|l3|) ). Moreover Sn(x) is not exactly ±1/2 even 
at x — when e is finite. Thus the above approximation 
becomes worse as the energy increases, making in general 
the magnitude of dI s /dV somewhat smaller than that of 
NfDQS. In particular the positive hump of dI s /dV at 
large V (near IOEd in this particular example) is smaller 
than the corresponding dip in Q near that energy. Hence 
at large voltages I s becomes negative as noted above. Q 

In conclusion, in this paper I studied the current- 
carrying density of states of a junction consisting of a 
dirty normal metal between two superconductors. I have 
also considered the dependence of the supercurrent be- 
tween the two superconductors on the applied voltages 
at the normal reservoirs of a cross geometry. 
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Appendix 

In this appendix I summarize some basic equations for 
easy reference, (c./., e.g. ||) The basic equation to be 
solved is the Usadel equation 

[er 3 ,g} + -d^(gd^) = (8) 

together with the normalization condition g 2 = —it 2 1 
governing the angular averaged matrix Green's function 
g which in turn has g R > A ' K , the retarded, advanced, 
and Keldysh martix Green's functions as its components. 
Here e is the energy. The pair-breaking mentioned in the 
text is simulated by e — > e + 17 where 7 > 0. [^2| 

g R can be parameterized as — i7r(cos#T3 — 
sinfteos^o^T! + sin9sm(f>a2T2) ■ g A can be related to g R 
by symmetry. The variables 9 and <j) obey the differential 
equations 

2i(e + i7)sin0 + D[d 2 x - sin6»cos6»(a0) 2 ] = (9) 
and 

d(sm 2 9d(t)) = (10) 

with the boundary conditions that they assume their 
equilibrium values at the reservoirs. For a normal reser- 
voir 9 = 0, while at a superconducting reservoir cos9 = 
-i(e + i-fj/V where V = y/A 2 - {e + ij) 2 . (7 -> 0+ if 
pair-breaking is not included). 

Q, related to the current-carrying density of states as 
discussed in the text, is given by 

Q = 2 Im[sm 2 9d(j)} (11) 

It is thus then position independent within any wire 
by eqn ([To]), a result which can also be directly ob- 
tained from the definition (^) for Q and by taking the 
appropriate trace of eq (|8|). Q obeys the symmetry 
Q(-e) = -Q(e). 

g K is expressed via the distribution function h as g R h— 
hg A where h can be chosen diagonal: h = /loTo + /13T3. 
The distribution functions obey the equations 
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d[Qho + (M 33 dh 3 + M 30 dh )] = 



(12) 



and the equation with <-> 3. These two equations ex- 
press respectively the conservation of particle and en- 
ergy at each individual energy (due to the absence of 
interactions). The (real) My coefficients are defined by 

The distribution functions at the reservoirs are given 
by their equilibrium values. At voltage V, ho(e) = 
[tanh^±§^ + tanh^f^]/2 and h 3 (e) = [tanh^f^ - 
tanh e "ty ]/2. Thus, at T — 0, when the voltage sweeps 
through the corresponding energy e = eV, the distribu- 
tion functions at y — —L y /2 change by Sho — — 1 and 
dh 3 = — 1. At the point where the voltage is —V ( 
y = L y /2 ), Sho = -1 and Sh 3 = 1. (see Fig. | ). 

The total number current density is given by 

N f D 



J N = 



de[Qh + (M 33 dh 3 + M 30 dh )} (13) 



The three terms represent respectively the contributions 
from occupation of current-carrying states, ordinary dif- 
fusion (with a modified diffusion coefficient) and an extra 
contribution due to broken particle-hole symmetry. 

If a potential barrier exists, there will be discontinuites 
of the parameters 9, 4> across the barrier. The appropri- 
ate boundary conditions are derived from |33| 

(2N f DS)gd^g = ±-[g( Xb _),g( Xb+ )} (14) 

where Rb is the resistance of the barrier at x b . 
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